
clear

bootidx = [1:80];

SaveRes = 'Results/full.mat';
load(SaveRes);

% Load Bootstrap estimates:
theta_boot_all = nan(14,length(bootidx));
p_uh_all       = nan(1, length(bootidx));
for b = bootidx
    bootname = ['Bootstrap/full_b', num2str(b), '.mat'];
    if exist(bootname) ~= 2; continue; end
    load(bootname)
    theta_boot_all(:,b) = theta_boot;
end

[theta_est, std(theta_boot_all,1,2)]



